ON ZEROS OF SELF-RECIPROCAL POLYNOMIALS 



MASATOSHI SUZUKI 



Abstract. We establish a necessary and sufficient condition for which all zeros of 
a self-reciprocal polynomial lie on the unit circle. Moreover, we relate the necessary 
and sufficient condition with a canonical system of linear differential equations (in the 
sense of de Branges). This relationship enable us to understand that the property of 
a self-reciprocal polynomial having only zeros on the unit circle is equivalent to the 
positive semidefiniteness of Hamiltonians of corresponding canonical systems. 



1. Introduction 

A nonzero polynomial P(x) = CQX n + c\x n ~ l + • • • + c n -\x + c n with real coefficients 
is called a self-reciprocal polynomial of degree n if Co 7^ and P(x) = x n P(l/x), or 
equivalently, cq 7^ and c& = c n _/% for every ^ k ^ n. The zeros of a self-reciprocal 
polynomial either lie on the unit circle, T = {z S C | \z\ = 1}, or are distributed sym- 
metrically with respect to T. Therefore, one of the basic problem is to find a "nice" 
condition of coefficients of a self-reciprocal polynomial for which all its zeros lie on T. 
In this paper, we study this problem for self-reciprocal polynomials of even degree. The 
restriction on the degree is not essential, because if P(x) is self-reciprocal and of odd 
degree, then there exists a self-reciprocal polynomial P(x) of even degree and an integer 
r ^ 1 such that P(x) = (x + l) r P(x). In contrast, the realness of coefficients are crucial 
to the results and methods of this paper. 

A nonzero polynomial P(x) = cox n + c\x n ~ x + • • • + c n -\x + c n (q € C) is called a 
self-inversive polynomial of degree n if cq 7^ and P(x) = fix n P(l/x) for some constant 
fx G C of modulus 1, or equivalently, cq 7^ and c& = nc n -k for every ^ k ^ n. If 
all coefficients of a self-inversive polynomial are real, then \jl must be ±1. Therefore, 
self-reciprocal polynomials are special cases of self-inversive polynomials. 

The Gauss-Lucas theorem |15] assert that the convex hull of the zeros of any non- 
constant complex polynomial contains the zeros of its derivative. Further, by a result of 
Schur |21j and Cohn [3], all zeros of a self-inversive polynomial P{x) lie on T if and only 
if all zeros of P'{x) lie inside or on T. Chen [1] found another necessary and sufficient 
condition for which all zeros of a self-inversive polynomial lie on T. It was rediscovered 
by Lahn-Smyth (See the book of Marden |16] and the survey paper of Milovanovic- 
Rassias |17| for several systematic treatments of zeros of polynomials.) The above results 
are often useful in order to examine whether all zeros of a given self-inversive (or self- 
reciprocal) polynomial lie on T, but a sufficient condition in terms of its coefficient is 
sometimes more convenient, because it may not be easy to check the above criterions for 
a given polynomial. 

Lakatos [S] presented a simple sufficient condition of coefficients of a self-reciprocal 
polynomial for which all its zeros lie on T. She proved that if P(x) G M.[x] is a self- 
reciprocal polynomial of degree n ^ 2 satisfying 

n-l 

c o| > ^2 \ck - co|, 
fc=l 
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then all zeros of P(x) lie on T (see Lakatos-Losonczi |10| . Kwon [6] and their references 
for further generalization of this kind of conditions) . Another simple sufficient condition 
in terms of coefficients was found by Chen [1] and Chinen [2], independently. They 
proved that if P{x) € R[x] has the form 

P(x) = (c x n + cix n ~ l + ■■■ + c k x n ~ k ) + (c k x k + c k _ x x k ~ x H + ax + c ), 

and Co > c\ > ■ ■ ■ > c k > (n ^ k), then all zeros of P{x) lie on T. As well as these 
results, known conditions of coefficients of general self-reciprocal polynomials for which 
all its zeros lie on T are usually sufficient conditions and they are not applied to the case 
that coefficients of middle terms are quite large comparing with coefficients of both ends. 

The first result of this paper is the construction of a set of 2g rational functions 
Rn(c) = Rg jn (c) G Q(c) (1 ^ n ^ 2g) of (g + 1) variables c = (co, • • • ,c g ) for every 
positive integer g such that all zeros of a self-reciprocal polynomial 

s-l 

P g (x) = J2 c k (x^~ k ^ + x k ) + c gX 3 (1.0.1) 

k=0 

with real coefficients c = (co, • • • , c g ) E ]R 9+1 (co 7^ 0) lie on T and simple if and only if 
R n (c) is a finite positive real number for every 1 ^ n ^ 2g (Theorem I2.6p . The rational 
functions R n (c) are constructed in Section [2.21 by using a linear system introduced in 
Section [2.11 The second result is a variant of the first result. We construct a set of 2g 
rational functions Rnic;^) = Rg t n(c; q 1 ^) € Q(c, q w ) (1 ^ n ^ 2g) of (g + 2) variables 
c = (co, • • • , c g ) and q u for every positive integer g such that all zeros of a self-reciprocal 
polynomial ([l.O.ip with real coefficients c = (cq, • • • , c 9 ) E (co 7^ 0) lie on T if and 

only if Rn{c\q u ) is a finite real number for every 1 ^ n ^ 2g and u; > 0, where q is an 
arbitrary fixed real number larger than one (Theorem I2.8p . Rational functions R n {c\q u ) 
are constructed in Section [2.31 by a way similar to R n (c). Moreover, it is proved that 

lim Rnic^q^ 1 ) = R n (c) 

as a rational function of c = (cq, • • • , c g ) (Theorem 12. 9p . 

Subsequently, we attempt to understand the above positivity conditions from a view- 
point of the theory of canonical systems of ordinary linear differential equations. After 
the review of the theory of canonical systems in Section [3j we construct two kinds of 
systems of ordinary linear differential equations attached to self-reciprocal polynomials 
in Section Bl (Theorem 14. 21 and I4.4p . It is proved that the system of Section |4~T1 (resp. Sec- 
tion I4.2p is a canonical systems if and only if all zeros of a corresponding self-reciprocal 
polynomial ([l.O.ip with real coefficients lie on T and simple (resp. lie on T) (Corollary 
WM (resp. Corollary IP]) ). 

We prove the result of Section 14.11 (Theorem I4.2p in Section [5] after a preparation of 
two basic tools. By using Theorem 14.21 and tools in Section [5j we prove the results of 
Section 12.21 in Section [6j Results of Section 14.21 are proved in Section [7] by a way similar 
to the proofs in Section [5] and [6l Finally, we comments on important remaining problems 
in Section [HJ 

This paper is written in self-contained fashion as much as possible. 

Acknowlegements This work was supported by KAKENHI (Grant-in-Aid for Young 
Scientists (B)) No. 21740004. 

2. Results I: 

TWO KINDS OF FUNCTIONS ATTACHED TO SELF-RECIPROCAL POLYNOMIALS 

In this section, we construct two kinds of rational functions R n (c) and Rnic]^) by 
using a linear system in order to state a necessary and sufficient conditions for all zeros 
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of a self-reciprocal polynomial to lie on T. Moreover, we clarify the relation between 
R n (c) and P n (c;<Z^)- Throughout this section, we fix a positive integer g. 



2.1. Linear System. Firstly, we define the (2k + 2) x (2k + 2) matrix P k (m k ) endowed 
with the parameter m k and the (2k + 2) x (2k + 4) matrix Q k for every nonnegative 
integer k as follows. For k = 0, 1, we define 



Pn := 





" 





1 



' 1 1 








1 -1 





" 1 













' l 





1 











Pi (mi) := 





1 








, Qi ■= 





1 




















1 














1 





-1 







1 





— mi 























For k ^ 2, we define P k (m k ) and Q k blockwisely as follows 





\v k + 







w k + 





Pk(m k ) := 





v k ~ 


j Qk '■ — 





w k 




04 


-m k ■ 0I k _ 




0fc,fc+2 


0fc,fc+2 



where 0I k := [0 k ,i h] , -m k 0I k = [0 k ,i -m k ■ I k ] , k [ is the k x I zero matrix, I k is 
the identity matrix of size k, 



v,: 



1 





























1 
















1 












1 










1 

























( 






















v 2 











1 





1 
























1 















1 





















" 







1 


















1 










1 












1 





( 











1 








1 






















1 


1 













1 























" 





1 




















-1 








1 














-1 














1 











-1 




















1 





-1 











1 





















" 







1 


















-1 










1 












-1 





( 











1 








-1 






















1 


-1 














X (k + 1) if A; is odd, 



fk + 2 



+ 1) if k is even, 



fc + 1 



x (k + 1) if k is odd, 



-±— j x (fc + 1) if k 



is even. 
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and matrices W k are denned by adding column vectors '(1 
end of matrices V^: 



0) to the right-side 



w: 



o 
i 








1 













1 







1 





1 











1 





k + 3 



x (k + 2) if k is odd, 



k + 2 



x (k + 2) if is even, 



WZ 



1 




















-1 





1 














-1 











1 











-1 

















1 





-1 











1 


















-1" 







1 












-1 













1 








-1 



















1 


-1 














k + 1 



(k + 2) if k is odd, 



k + 2 



x (k + 2) if k is even. 



Lemma 2.1 For k ^ 1, we have 



det P k (m k ) 



e 2j +i 2 J m^ x ifk = 2j + l, 



£2j 2 3 m 



2j 



where 



£ 2j+l 



+1 j = 2,3 mod 4 
-1 j = 0, 1 mod 4 



£ 2j 



ifk = 2j, 



+1 j = 0,1 mod 4 
-1 j = 2,3 mod 4 



In particular, P k (m k ) is invertible if and only if m k ^ 0. 

Proof. This is trivial for k = 1. Suppose that = 2j + 1 ^3 and write P k (m k ) as 
(«l • • • U2fc+2) by hs column vectors uj. At first, we make I k+2 a t the left-upper corner 
by exchanging columns V( k+5 y 2 , ■■■ , Vk+i and v k+2 , ■■■ , f (3fc+3)/2 so that 

detP fc (m fc ) = det(ui • • • V( k+3)/2 v k+2 ■ ■ ■ v^ k+3)/2 V( k+ 5)/2 • ■ • v k+1 «( 3 fc+5)/2 • • • «2fc+2)- 
Then, by eliminating every 1 and — m k under I k + 2 of the left-upper corner, we have 



det P k (m k ) = det 



h+2 * 

0fc,fc+2 Z k 



Here Z k is the k x k matrix 



Z k 



Z k ,i z k2 
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for which Zk i is the j x j anti-diagonal matrix with —1 on the anti-diagonal line and 



U 


u 




U 


-m k 






n 
U 




U 








-m k 


-m fc 









-m k 
















-2m k 
















-2m k 




















—2m k 



j x U + 1) 



(j + 1) X (j + 1) 



The above formula of det P k {m k ) implies the desired result. The case of even k is proved 
by a way similar to the case of odd k. □ 

Let z = (xo, ■ ■ ■ , X2g,yo, • • ■ > U2g) be {Ag + 2) indeterminate elements. We define the 
column vector v g (0) = v g (z; 0) of length (4g + 2) by 

V g (0)= t (x ■■■ X 2g V0 ■■■ V2g) , 

where v means the transpose of a row vector v, and define column vectors v g (n) = 
v g(z; n) (1 ^ n ^ 2g) of length (4g — 2n + 2) inductively as follows: 

tfc(ra-l)[l]+t; fl (n-l)[20-n + 2] 



m 2g-n(l) - 



v g (n - l)[2g - n + 3] - u 9 (rt - l)[4g - 2n + 4] ' 



U 9 (n) := P2 5 -n("l2 9 -nU)) 1 Q2g-nV g {n - 1), 



(2.1.1) 
(2.1.2) 



where Po( m o) := -fb an( l v \j] means the jth component of a column vector v. By the 
following lemma, we confirm that m2g- n (z) and v g (n) are well-defined as a rational 
function and a vector valued rational function of z = (xo, • • • , X2 g ,yo, ■ ■ ■ , V2g) f° r every 
1 ^ n ^ 2g, respectively. 

Lemma 2.2 Let K = Q(xq, • • • , X2 g ,yo, ■ ■ ■ , V2g)- For every 1 ^ n ^ 2g, m2g- n (z) is a 
nonzero element of K and hence P2g~ n ( m 2g~n(z)) € GL^g^ 2 n+2(K). 

Proof. By Lemma 12.11 it is sufficient to show that there exists a numerical vector c G 
R 9+2 such that |m2 9 - n (c)| < oo and ra2 9 - n (c) 7^ for every 1 ^ n ^ 2y. Existence of 
such numerical vector c is guaranteed by the necessity part of Theorem 12.41 below, since 
it is clear that there exists a self-reciprocal polynomial of degree 2g with real coefficients 
having only simple zeros on T for every positive integer g. □ 

Here we mention that the vector v g (n) of (|2.1.2p can be defined from v g (n — 1) by a 
slightly different way according to the following lemma. 

Lemma 2.3 For every 1 ^ n ^ 2g, we have 

Vg(n-l)[i\+v g (n-l)[2g-n + 2] 



v g (n - l)[2g - n + 3] - v g (n - l)[Ag - 2n + 4] 

{P2g-n{m2g-n)~ 1 Q2g-nVg{n - 1))[1] 



(P2g-n{m2g-n) 1 Q2g-nV g (n - l))[2g -71 + 2]' 

that is, the right-hand side is independent of the indeterminate element m2 g - n - 
Proof. Define (k + 1) x (k + 2) matrices M kfl , M fe)2 , M fc;3 , M fc , 4 by 

"l 10 

110 1 [Mi,! Afi )2 ] 1 

1-1' Afi! 3 Mt! 4 1 

1/mi 



"M ,i 


M , 2 " 




M ,3 


M ,4_ 
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and 



k,l 



2 
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o 
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n 
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o 




n 
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o 


n 
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n 
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n 
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n 




n 
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n 
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2 
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1 




















1 


1 
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1 











1 


















1 






if k ^ 3 is odd, 



if k ^ 2 is even, 



M K 



2 



2 
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-1 


















1 





-1 



































if 










-1 





1 


















-1 
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-1 
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1 
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-1 
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-1 






















-1 
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-1 
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1 






if k ^ 3 is odd, 



if k ^ 2 is even, 
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if k > 3 is odd, 



if k ^ 2 is even, 
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1 



if k ^ 3 is odd, 



if k ^ 2 is even. 



Then, we find that 



Pk{m k ) 1 Q k 



M k>1 



M K2 



(2k + 2) x (2k + 4) 
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for every ^ k ^ 2g. This formula of Pk(mk) l Qk shows that 

(P2g-n(m 2 g-n)~ 1 Q2g-nVg(n - 1))[1] 

= v g (n - 1)[1] + ^(n -l)[2g-n + 2], 

(P2g-n{'m2g-n)~ 1 Q2g-nVg{n - l))[2g -71 + 2] 

= v s (n - 1)[2# - n + 3] - v g (n - l)[4g - 2n + 4] 
for every 1 ^ n ^ 2g. These equalities imply the lemma. □ 

By Lemma [2. 31 we can define v g (n) by taking v g (n) = P 2 g-n{fn2g-n)~ l Q2g-n v g {n—l) 
for v g (n — 1) and then substituting the value v g (n)[l]/v g (n)[2g — n + 2] into v g (n). 
Anyway, staring from the initial vector v g (0), 2g vectors 

v g (l), v g (2),--- ,v g (2g-l), v g (2g) 

are defined by using Pjfc(mjt) and Q^. They have entries in Q(xq, ■ ■ ■ , x g , yo, ■ ■ ■ ,y g ), in 
other words, their entires are functions of the initial vector v g (0). 

2.2. The first result. Let c = (co, • • • , c g ) be (g + 1) indeterminate elements and let q 
be a real variable. We take the column vector 



w ff (0) 



a 9 (0) 

W(o) 



,(o) 



CO 
Cl 

=9-1 

C 9 
=9-1 

Cl 
CO . 



b s (0) 



co log(g 9 ) 
Cl log(<^-i) 

c 9 _ilogg 


-c 9 _ilogg 

-diog^- 1 ; 

-co log(g 9 ) 



(2.2.1] 



of length (4g + 2) as the initial vector of the system consisting of (|2.1.ip and (|2. 1 .2[) . 
Then 2g column vectors 

v g (l), v g (2),--- ,v g (2g-l), v g (2g) 

are defined with enties in Q(ci,--- ,c g ,logq) and m 2g - n (c ; log q) € Q(ci, 
are defined by 

v g (n-l)[l]+v g (n-l)[2g-n + 2} 



,Cg, lOg ( 



"l2g-n(c;logg) : 



v g (n - l)[2g - n + 3] - - 1)[4# - 2?i + 4] 

(P 2g -n(m2 g -n)~ 1 Q2g-nVg(n - !))[!] 
(P2g-n(m2g-n)~ 1 Q2g-nVg(n ~ l))[2g - U + 2] 

v 3 (n)[2 5 -n + 2] 



(2.2.2) 



for 1 ^ n ^ 2g according to (|2.1.ip and (|2.1.2p . where the second equality follows from 
Lemma 12.31 and the third equality is definition (|2.1.2p . 
Now the first result is stated as follows. 



Theorem 2.4 Let g 1 and q > 1. Let P g (x) be a self-reciprocal polynomial (jl.O.ip 
of degree 2g with real coefficients c = (co, ■ ■ ■ , c g ). Define 2g numbers m 2g - n (c',logq) 
(1 2g) by substituting the numerical vector c into (|2.2.2|) . Then all zeros of P g {x) 

lie on T and simple if and only if 

m 2g -n{c]\ogq) > and m 2 g- n (c ; log q)~ l > 

for every 1 $C n ^ 2g. 

Proof. This is proved in Section [6] □ 
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Remark The necessity part of Theorem 12.41 shows that v g (n) and ra2 9 - n (c, ;logg) are 
well-defined for the special initial vector (|2.2.ip as well as in the proof of Lemma 12.21 



Remark The criterion of Theorem 12.41 can be written as rri2g-n (c ; log q) > and 
m 2g-n{c ; log q)~ l ^ 0. The reason why we state it as in Theorem 12.41 is for the conve- 
nience of later subjects. See (|4.1.6p . (|4.1.7p and Corollary 14.31 below. It is similar about 
Theorem EH EZ3 El below. 

Remark The strict inequalities mi g — n {c ; log q) > and m2 5 _ n (c ; log g)" 1 > are 
essential. In fact, the self-reciprocal polynomial P2(x) = 4(x 4 + 1) — 16(x 3 + x) + 23x 2 = 
(2x 2 — 3x + 2)(22; 2 — 5x + 2) has a zero outside T but 1713 = (2 log g) -1 > 0, m 2 = m\ = 0, 
mo = (126 log g) -1 > for any q > 1. 

The criterion of Theorem 1 2. 41 is independent of the choice of q > 1. In order to clarify 
such thing, we modify the statement of Theorem 12.41 as follows. 
In addition to definition (|2.2.2p . we set the convention 

m 2g (c;logq) := — — — . (2.2.3) 
glogq 

We define R n (c) = R g , n {c) for < n < 2g by 

Rq(c) := 1 



and 



Rn(c) :-- 



TT m 2g~(2j+l){c;logq) 

7=0 m 2<?~2i(c;logg) 

n m2g - (2j+2) i g; ;° gg ; if„= 2 j +2 >2 

v f= m 29-{2j+i){c;\°gq) 



if n = 2J + 1 > 1, 

(2.2.4) 



By the definition, we have 

m2 9 - n (c ; log q) = — — — R n -i{c)R n (c) (2.2.5) 
5 log 9 

for every 1 ^ n ^ 2g. 

Lemma 2.5 We /iawe R n (c) G Q(co, • • • , c ff ) /or every ^ n ^ 2g. 

Proof. Put F = Q(c , • • • ,c 9 ). At first, we show that v s (ra)[fc] G F for 1 ^ k < 2g-n + l 
and (log <7)~ 1 u 9 (n)[A;] € F for 2g — n + 2 ^ k ^ 4o — 2n + 2. It is clear for v g (0) by 
definition (|2.2.ip . Assume that it holds for v g (n — 1). Then, by definition (|2.2.2p . we have 
w-2 5 -n(c;logg) = (logg)~V2 9 -n(c) for some fJ, 2g -n(c) G F. By applying the formula 
of Pk{mk)~ l Qk m the proof of Lemma 12.31 to k = 2g — n, we obtain v g (n)[k] G F for 
l^/c^2g-n+l and (logg)" 1 t; !? (n)[A;] G F for 2g - n + 2 < k < 4g - 2n + 2. Hence 
(log q)m2g- n {c ; log g) G F for every 1 ^ n ^ 2g by induction. 

By definition (|2.2.4p . R n (c) has the same number of m^g-jic. ; log g)'s in the denomi- 
nator and the numerator. Hence R n (c) G F for every ^ n ^ 2n. □ 

Let c G and g > 1. If R n (c) > for every 1 ^ n ^ 2g, then ?n2 9 _ n (c ; log g) > 

for every 1 ^ n ^ 2<? by (|2.2.5p . Conversely, if m2 9 _„(c ; log g) > for every 1 ^ n ^ 2g, 
then Fi n _i(c) and R n (c) must have the same sign for every 1 ^ n ^ 2g by ()2.2.5p . 
However, Po(c) = 1 > by the definition. Hence R n {c) > for every 1 ^ n ^ 2g. As a 
consequence, Theorem 12. 41 is equivalent to the following statement. 

Theorem 2.6 Let P g {x) be a self-reciprocal polynomial (jl.O.ip of degree 2g > wii/i 
reaZ coefficients c = (cq, • • • , c s ). Define 2g numbers R n (c) by substituting c into (|2.2.4p . 
TTien a// zeros of P g (x) lie on T and simple if and only if 

Rn(c) > and R n {c)~ l > 

for every 1 ^ n ^ 2g. 
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For small positive integer g, we can calculate R n (c) by hand according to the defini- 
tion. Because i?i(c) = 1 by the definition, interesting values are -R 2 (c), ■ ■ ■ , R2g(c)- F° r 
example, we have the following small table: 

ID/ > 2c + c l 

•5 = 1) #2(co,ci) 



2c - ci 

g = 2, R n = Rn(cQ,ci,c 2 ) (2 < n < 4), 



i?2 



4c + ci 



8cg 



2cf + 4c c 2 



4c - ci ' 8cq + c\- 4c c 2 

5 = 3, R n = i? n (co,ci,c 2 ,c 3 ) (2 < re < 6), 
6cq + ci n 18cq - 3cf + 6c c 2 



R i 



2c + 2ci + c 2 
2c - 2ci + c 2 



i?2 



6c - ci 

,3 



18 C 2 + 2cj - 6c c 2 : 

.3 



R 4 



36cg + 6c^ci - c cf + 4cf - 14c cic 2 + cfc 2 - 4c c| + 18c^c 3 + 3c cic 3 
36qj - 6cqCi - c cf - 4cf + 14c cic 2 + cfc 2 - 4c c| - 18cqC 3 + 3c cic 3 ' 
R 5 = (1084 ~ 21c o c i - 12c i + 108cgc 2 + 42c c?c 2 

— 12cqC 2 + 3cfc 2 — 12coc 2 — 54cqCic 3 — 6cfc 3 + 30coCic 2 c 3 — 27cqc|)/ 
(108cg + 9c§c? + 8cj ~ 108c|]c 2 - 42c c?c 2 

+ 36cqc| + cic 2 - 4c c 2 + 54cqCic 3 - 4cfc 3 + 18c cic 2 c 3 - 27cqC 3 ), 
2c + 2ci + 2c 2 + c 3 

-Kg — . 

2c - 2ci + 2c 2 - c 3 
See also Section [8] (1) for another formula of R n (c). 

2.3. The second result. In order to deal with the case that P g (x) may have a multiple 
zero on T, we consider the following variant of m 2s _ n (c ; log q) and R n (c). 



In stead of the vector (|2.2.ip . we take the column vector 



Vg(0) 



a <7,w(°) 



(0) 



(0) 



co q 



ci q 



Cg-iq 
Cg-\ q- 



-gu 



(2.3.1) 



of length (4g + 2) as the initial vector of the system consisting of (|2.1.ip and (|2.1.2p . 
Then 2g vectors 

^(1), v g (2),--- ,v g (2g-l), v 9 (2g) 
are rational functions of c = (cq, • • • ,c g ) and q u over Q. We define 

v g (n-l)[l]+v g (n-l)[2g-n + 2] 



m 2 g- n (c;q UJ ) : 
for 1 ^ re ^ 2g and 



v g (n - l)[2g - n + 3] - u 9 (n - 1)[4# - 2n + 4] 



m2 g {c;q ul ) :-- 



q 9" + q -gu 



qgu — q-gu 

as well as (f2~2~2l and (pEOD . Further, we define i2 n (c; for < re < 2# by 

i? (c;O :=1 



(2.3.2) 



(2.3.3) 



ON ZEROS OF SELF-RECIPROCAL POLYNOMIALS 



11 



and 



R n (c;q u 



^rn 2g ^ +1) (c-, q p . {n = 2J + 
TT m23 - (2j+2) |- ;g 1 1 fn = 2J + 2^2. 



(2.3.4) 



Then, we obtain the following second results. 

Theorem 2.7 Let q > 1. Let P g (x) be a self-reciprocal polynomial (|1.U.1|) of degree 2g > 
with real coefficients c = (co, • • • , c g ). Define 2g numbers m 2g _ n (c; q u ) by substituting 
c into (|2.3.2p . Then, all zeros of P g (x) lie on T (allowing multiple zeros) if and only if 

m 2g - n (c;q w ) > and m 2g _ n (c;q UJ )~ 1 > 
for every 1 ^ n ^ 2g and uj > 0. 

Theorem 2.8 Let q > 1. Let P g (x) be a self-reciprocal polynomial (jl.O.ip of degree 
2g > with real coefficients c = (co, • • • , c fl ). Define 2g numbers R n {c ; q^) by substituting 
c into (|2,3,4p . Then, all zeros of P g (x) lie on T (allowing multiple zeros) if and only if 

R n (c;q u) )>0 and R n (c;q u )~ l > V 

for every 1 n 2g and u > 0. 

The rational functions of (|2.2.4p and (|2.3.4|) have the following simple relation. 

Theorem 2.9 Let i?„(c;<f ) and Rn(c;q u ) be in (f2~2~4")) and (|2.3.4|) . respectively. Then 

lim R n (c;(f)=Rn{c) 

as a rational function of c = (co, • • • , c s ) over Q. Suppose that all zeros of a self-reciprocal 
polynomial (jl.O.ip with real coefficients c = (co, • • • , c ff ) lie on T and simple. Then we 
have 



R n (c;q UJ ) = Rn(c) + 0(logq UJ ) as q" -> 1" 



and 



m 2g - n (c;q u ) = -^—(Rn-i{c)Rn{c) + 0(logq u )) +0(logg w ) as <f -»■ 1+, 
glogg^ V / 

= -(m 2s - n (c; log g) +O(w/0)) +0(w log g) as w -> + , 

where implied constants depend only on c. 

Theorem 12.71 12. 8| and 12.91 are proved in Section [7] together with Theorem 14.41 and 
Corollary 14.51 below. 

3. Canonical systems 

The positivity of values m 2ff _ n (c; logo), m 2g - n (c; q u ), R n (c), R n (c; q w ) for a numeri- 
cal vector c attached to a self-reciprocal polynomial of degree 2g having only (simple) ze- 
ros on T can be understood from the viewpoint of canonical systems of linear differential 
equations (in the sense of de Branges). In fact, ideas of constructions of m 2g - n (c ; log q) 
and m 2g - n (c;q UJ ) in Section [2] are coming from the theory of canonical systems. 

In this section, we review the theory of canonical systems of linear differential equations 
and the theory of entire functions of the Hermite-Biehler class according to de Branges 
[1], Dym [5], Levin |13| . Remling |18j . and Lagarias 0(8], in order to understand Theorem 



2.41 and 12.71 in terms of these theories. We often use the notation F^(z) = F{z) for an 
entire function F(z). An entire function F(z) is called real or a real entire function if 
F(M) C K, or equivalently, F*(z) = F(z). We denote by lim x + and lim ^ - the 
right-sided limit and the left-sided limit at x = xq, respectively. 
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Definition 3.1 Let H(a) be a 2 x 2 matrix-valued function defined almost everywhere 
on a finite interval I = [ctx,«o) ( — 00 < a i < a < °°)- ^4 family of linear differential 
equations on I of the form 



~A(a, z)~ 
B(a,z) 


= z 



1 


-1 



H(a) 


~A(a,z)~ 
B(a,z) 


, lim 

a— »a^~ 


"A(a, z) 
B(a,z) 




"1* 




parametrized by z € C is called a (two-dimensional) canonical system if 

(HI) -ff(a) is a positive semidefinite symmetric matrix for almost every a £ I, 
(H2) i-f(a) ^ on any open subset of I with positive Lebesgue measure, 
(H3) -ff(a) is locally integrable on I. 

For a canonical system, the matrix-valued function H(a) is called its Hamiltonian. 

Usually, canonical systems are defined by using the additive derivative d/dt for t = 
log(ao/o)- However, we use the multiplicative derivative —a(d/da) as above for the 
convenience of descriptions of results in the paper. 

Canonical systems are closely related to entire functions of the Hermite-Biehler class. 

Definition 3.2 (Hermite-Biehler class) An entire function E{z) is said to be a function 
of class HB if it satisfies the condition 

\E\z)\ < \E(z)\ for every lmz>0 (3.0.5) 

and has no real zeros. On the other hand, an entire function E(z) is said to be a function 
of class HB if it satisfies the condition 

\E^(z)\ ^ |^(^)| for every Imz > 0. 
and has no zeros in the upper half-plane Imz > 0. 

Remark This definition of class HB is equivalent to the definition of Levin |13[ §1 
of Chap. VII] if we replace the word "the upper half-plane" by "the lower half-plane", 
because (|3.0.5p implies that E(z) has no zeros in the upper half-plane Imz > 0. We 
adopted this definition for the convenience of using of the theory of canonical systems 
via the theory of de Branges spaces. 

The following result for a function of class HB is used often in the later sections. 

Proposition 3.3 Let E{z) be an entire function of finite order. Put 

A{z) := \(E(z) + Efl(z)), B(z) := % -{E{z) - E\z)). 

Then E(z) is a function of class HB if and only if E(z) has no zeros in the upper half- 
plane Imz > 0, (real) entire functions A(z) and B(z) have only simple real zero, and 
zeros of A(z) and B{z) interlace. On the other hand, the function E(z) is a function of 
class HB if and only if E{z) is a product of a real entire function Eq(z) having only real 
zeros and a function E\(z) of class HB. 

Proof. See Levin |13| Chap. VII, Theorem 3, Theorem 5, the latter half of p. 313]. □ 

There are two important results of de Branges that relate a canonical system with 
an entire function of class HB. Roughly, if (A(a, z), B(a, z)) is a solution of a canonical 
system, then E(a, z) := A(a, z) — iB(a, z) is a function of class HB for every [a±, ao) (see 
[3] Theorem 41], and also Section 2] for details). Conversely, if E{z) is a function 
of class HB normalized as E(0) = 1, then there exists a canonical system on some 
interval [ai,ao) such that E(z) = A(ai,z) — iB(a\,z) for the solution (A(a, z), B(a, z)) 
of the canonical system (see [H Theorem 40], and also |18[ Theorem 7.3], pp. 70-71] for 
more details). Therefore, a function of class HB and a Hamiltonian of a canonical system 
correspond each other. However, in general, an explicit construction of a Hamiltonian is 
quite difficult when we start from a function of class HB, because it is a kind of inverse 
scattering problem. 
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4. Results II: 

Differential equations attached to self-reciprocal polynomials 

In this section, we construct two kinds of systems of ordinary linear differential equa- 
tions attached to a self-reciprocal polynomial so that the first one (resp. the second one) 
is a (two dimensional) canonical system if and only if all zeros of the polynomial lie on 
T and simple (resp. lie on T). 

4.1. Systems of the first kind. Fix a real number q > 1 arbitrary. For a self-reciprocal 
polynomial P g (x) of (jl.O.ip with real coefficients c = (cq, ■ ■ ■ ,c g ), we define 

9-1 

A q {z) := q~ 9lz P g (q lz ) = £>(« to - fc) " + q~^ iz ) + c g . (4.1.1) 

fc=0 

By the definition, all zeros of P g (x) lie on T if and only if A q (z) has only real zeros. The 
self-reciprocal condition P g {x) = x 2g P g (l/x) implies the functional equation A q (z) = 
A q (—z) and the realness of coefficients of P g (x) implies A q (z) = A q (z). Hence A q (z) is 
an even real entire function of exponential type. Further, we define 

B q {z) := —Mz) (4.1.2) 

and 

E q (z):=A q (z)-iB q (z). (4.1.3) 
Then B q (z) is a real entire function and 

g-i 

B q (z) = -t(Iog?) ]>> - k)c k (q^ z - q-(9-^y 

k=0 

Moreover, we have 

El(z) = A q (z)+iB q (z) (AAA) 

and 

M*) = \{E q {z) + E\{z)), B q (z) = l -(E q (z) - E\{z)). 

Lemma 4.1 Let E q (z), A q (z), B q (z) be as above. Then 

(1) E q (z) satisfies condition f|3.0.5|) if and only if A q (z) has only real zeros. 

(2) E q (z) is a function of class HB if and only if A q (z) has only simple real zeros. 

By this lemma, if A q (z) has only real zeros, then E q (z) is a function of class HB at 
least and B q (z) has only real zeros, but A q (z) and B q (z) may have a common (real) zero. 

Proof. (1) Assume that E q (z) satisfies (|3,0.5p . Then it implies that A q (z) ^ for 
Imz > 0, Further, A q (z) ^ for Imz < by the functional equation A q (z) = A q (—z). 
Hence all zeros of A q (z) lie on the real line. Conversely, assume that all zeros of A q (z) 
are real. Then A q (z) has the factorization 

A q {z) = C\im JT (C, pGR), (4.1.5) 

_R->oo - L - L V pi 

because A q (z) is real, even and of exponential type. Therefore, 

Hence, for Imz > 0, 

\E q (z)\ = \A q (z)\ 



1 + i 



AJz) 



> \M*) 



Mz) 



\E\{z)\. 
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(2) Suppose that E q (z) is a function of class HB, that is, E q {z) satisfies (|3.0,5p and has 
no real zeros. Then, A q {z) has only real zeros by (1). If A q (z) has a multiple real zero, 
then A q {z) and B q (z) = —A' (z) have a common real zero. Thus E q (z) = A q (z) —iB q {z) 
has a real zero. It is a contradiction. Hence A q (z) has only simple real zeros. The 
converse assertion follows from (1) and definition (|4.1.2p . □ 

As mentioned in Section [3J functions of class HB and canonical systems correspond 
each other. Therefore, by Lemma 14.14 there must exists a canonical system attached to 
E q {z) if and only if all zeros of P g (x) lie on T and simple. Now we construct a system 
of linear differential equation so that it is a canonical system if and only if all zeros of 
P g (x) lie on T and simple. 

Let c = (co, • • • , c g ) be 2g indeterminate elements and let q > 1 be a real number. 
By using m2 g - n (c ; log q) of (|2.2.2p . we define the Q(co, • 
a€ [l,q 9 ) by 

n — 1 

m q {a) = m 2 g- n (cAogq) if q~ ^ a < 
and the 2x2 matrix valued function H q (a) of a € [1, q 9 ) by 



Cg)- valued function m q (a) of 



(4.1.6) 



H q (a) 











In addition, we define functions A q (a, z) and B q (a,z) of (a, z) £ [1,9 s ) x C by 



T n (a,z)v g (n) if q ? ^a<q? 



A q (a,z) 


1 


"l " 


B q (a,z)_ 


~ 2 


-i 



c g {a,z) 




where T n (a, z) is the 2 x (4g — 2n + 2) matrix valued function 

c g -i(a,z) ■■■ c^ g+n (a,z) 

••• s g (a,z) s g -i(a,z) ■ 

with 

c k (a, z) := {{q k /aT + {q k /a)-% s k (a, z) := ((q k / a y z - 
Then the result of this section is stated as follows. 







(4.1.7) 



(4.1.8) 



(4.1.9) 



•• s_ 3 _|_ n (a, z) 
(q k /a)- iz ). (4.1.10) 



Theorem 4.2 Let c = (co,-- - ,c g ) be real coefficients of a self -reciprocal polynomial 
Pg{x) of degree 2g in (jl.O.ip and let q > 1 be a real number. Define E q {z), A q {z), 
B q {z) by ([4. 1.31) . (|4.1.1|) . (|4.1.2p . respectively. Define H q (a), A q (a,z) and B q (a,z) by 
substituting the numerical vector c into (j4 . 1 . 7|) . (|4.1.8p . respectively. Define 2g numbers 
r n2g-n(c',logq) (1 ^ n ^ 2g) by substituting c into (|2,2.2p . Then we have 

(1) m 2 g-i(c;loggr) >0 and m 2g -i(c ; log g)" 1 > 0, 

(2) ^4^(0,2;) and B q (a, z) are defined on [l,q 1//2 ) and 

A q (l,z)=A q (z), B q (l,z) = B q (z). 

Let 1 no ^5 2g. Suppose that m2g~ n (c;logq) 7^ and m2 g - n (c ; log q)" 1 7^ /or euery 
1 ^ re ^ no- TTien we Ziaue 

(3) A q (a,z) and B q {a,z) are defined and continuous on [l,q n °l 2 ) with respect to a, 

(4) the left-sided limit lim a _^/ g „ /2\_ (A q (a, z), B q (a, z)) defines entire functions of z, 

(5) A q {a,z) andB q (a,z) are differentiate functions on (g( n-1 )' 2 , q n l 2 ) with respect 
to a for every 1 ^ n ^ hq, 

(6) H q {a), A q (a,z) and B q (a, z) satisfy the system 



d_ 

da 



A q (a,z) 
B q (a,z) 



H q {a) 



A q (a,z) 
B q (a,z) 



(z 6 C) 



(4.1.11; 



for a G [1, q n ° /2), 

Suppose that A q (a, z) and B q (a,z) are defined on [q 9 ~ 1 ^ 2 ,q 9 ), that is, substitution of c 
into v g (2g) defines a numerical column vector of length 2. Then, we have 
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(7) lim 

a-Kg»)- 



A q (a,z) 
B q (a,z) 



A(0) 

Proof. This is proved in Section [5j 



E q (0) 



□ 



Corollary 4.3 The 2x2 matrix valued function H q (a) of (|4.1.7p attached to real coef- 
ficients c of P g (x) defines a canonical system on [l,q 9 ) if and only if all zeros of P g (x) 
lie on T and simple. If H q (a) defines a canonical system on [l,q 9 ), then the pair of 
functions (A q (a, z)/E q (0), B q (a, z) / E q (0)) is a solution of the canonical system. 

Proof. If all zeros of P g (x) lie on T and simple, then A q (z) has only simple real zeros. In 
particular, E q (0) = A q (0) ^ 0, since A q (z) is even and B q {z) is odd. Hence, by Theorem 
12.41 and Theorem 14.21 H q (a) defines a canonical system on [l,q 9 ) and (A q (a, z)/E q (0), 
B q (a, z)/E q (0)) is its solution. Conversely, if H q {a) defines a canonical system, then 
m 2g-n{c;logq) > and m2 9 _ n (c ; log q)~ l > for every 1 ^ n ^ 2g by Definition 13.11 

□ 



Hence, all zeros of P g (x) lie on T and simple by Theorem |2~1 

4.2. Systems of the second kind. By using the function of (|4.1,ip . we define 

E qju (z) := A q (z + iu), (4.2.1) 

A q>ul (z) := \{E q ^(z) + 4, w (*)), B q ^(z) := l -{E q ^(z) - E^z)). 
Then A q ^{z) and B q ^{z) are real entire functions satisfying 

E q ,uj(z) = A qi0J (z) - iB q ^(z) 



and 

Therefore, we obtain 

A q ,Lu(z) 



E*Jz) = A q (z - iu). 



- (A q (z + iu) + A q (z - iu)) 



k=l 



kuj _|_ —ku) \ I „kiz 



+ 



-kiz 



+ C 



y- 



B q ,U z ) = - (A q (z + iu) - A q (z - iu)) 



2 

k=l 



ku) \ ( „kiz 



-kiz 



Further, by using v g (n) and rn2 g - n {c\q u ) of Section f2. 31 we define the function m q>UJ (a) 
of ae [l,q 9 ) by 

v g (n)[l] 



m 2g ^n(c;q u ') :-- 



v g (n)[2g - n + 2} 



if q 2 ^ a < q 2 



and define the 2x2 matrix valued function H qiU) (a) of a G [1, q 9 ) by 
In addition, we define functions A qtU) {a, z) and B q ^(a, z) of a € [1, q 9 ) by 



m q ^{a) 1 
m q ^(a) 



Aq t cj{a, z) 


1 


"1 " 


_B q!U] (a,z)_ 


~ 2 


-i 



T n {a,z)v g (n) if q 2 ^ a < q^ 



as well as ()4.1.8p . where T n (a, z) is the 2 x (4g — 2n + 2) matrix valued function in (|4.1.9p . 
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Theorem 4.4 All assertions of Theorem \4-%\ hold if we replace E q (z), A q (z), B q (z), 
H g (a), A q (a,z), B q (a,z) by E q ^(z), A q ^(z), B q ^(z), H q ^(a), A q ^(a,z), B q ^(a,z), 
respectively. 

Corollary 4.5 The 2x2 matrix valued function H q ^ w {a) defines a canonical system on 
[l,q 9 ) for every u > if and only if all zeros of P g (x) lie on T. If H q uJ (a) defines a 
canonical system, then {A q ^{a, z)/E q ^(0), B q ^{a, z) / E q ^{Q)) is its solution. 

We prove Theorem 14,41 and Corollary 14,51 in Section [71 



5. Proof of Theorem 14.21 
Proof of (1) and (2). By definition fl2~2~2]) . 

VO)[l]+^(0)[2<7 + l] 1 



m 2g -i 



v g (0)[2g + 2] - v g (0)[4g + 2] glogq 



> 0, 



(5.0.2) 



since q > 1. This implies (1). By definition (|2.1.2p . v g {l) = P 2g -\{m 2g ^i) 1 <?2g-i '^g(O). 
Therefore, we have 



%(1) 

Vi). 



a g (l) 



2c 

ci(l + m 2 g_ilog(^~ 1 )) 
c 2 (l + m 2g -i\og{q9- 2 )) 

c 3 _i(l + m 2 g-ilogg)) 
Cy-i{\ - m 2g -i\ogq)) 

c 2 (l - m 2g -i\og{q 9 - 2 )) 
ci(l - m 29 „ilog(g9~ 1 )). 



By substituting (|5.0.2|) into this formula of v g (l), we obtain 



b ff (l) 



2c log(^) 

c l( m 2g-l + lo g(9 9_1 )) 
c 2 {m' 2g \ l + \og{q3- 2 )) 

C ff-l( m 2g-l + lo g<?) 

-1 

C 9 m 2g-1 

Cg-l^p-l - lo S<?) 

c 2(^2 ff -l -log(^~ 2 )) 
C l( m 2g-1 ~ \og(q 9 ~ 1 ))_ 



%(1) 
b,(l) 



M 1 ) 



2c 
2g-l 


2.9-2 



ci 

C2 



2±i C 1 
9 C 9- 1 

9-1 



9 Cl 



Mi) 



2#c log g 
(25 - l)ci logg 
(2 5 -2)c 2 logg 

(§ + l)c g -ilogg 

0c 9 log g 
(5 - l)c g -i log g 

2c 2 log g 
ci log q 



(5.0.3) 



By definition (j4~l~8j) . 



~A q (l,z)' 


1 


"1 " 


B q (l,z)_ 


~ 2 


-i 



Cg(l,z) Cg-l(l,Z) 





-<?+l(M) 






Sg(l,Z) Sg-l{l,z) 





s- g+1 (l,z) 



Vg(l). 



ON ZEROS OF SELF-RECIPROCAL POLYNOMIALS 



17 



On the other hand, it is easy to see that C- k (l, z) = c k (l,z), S- k (l,z) = — s k (l,z), 
co(l, z) = 2, so(l, z) = by definition ()4.1,10p . Hence we obtain 

A q (l, z) = y~] c k c g - k (l, z) +c g = A g (z), 

k=0 

9-1 

Bq(l, z) = -z(log q) ^(g - k)c k s g - k (l, z) = B q (z) 

k=0 



and complete the proof. 



□ 



Proof of (3) and (4). By the assumption and definitions (|2,1.2p and (|2,2,2p . numerical 
vectors v g (l),--- ,v g {no) are well-defined. Thus, functions A„(a, z) and B q (a,z) are 
defined on [l,q n °/ 2 ) and (4) is trivial by definition (|4.1.8p . Therefore, it is sufficient to 
prove that A q (a, z) and B q (a, z) are continuous at a = g( n_1 )/ 2 for every 1 ^ n ^ no — 1. 
By definition (|4.1.8p . the continuity of A q (a, z) and B q (a, z) at a = ^("-i)/ 2 is equivalent 
to the equality 

Tniq^-W 2 , z)v g (n) = T^q^' 2 , z)v g {n - 1), 
and this is also equivalent to 

T n { q ^l\z)P 2 g^ n {m 2 g^ n r 1 Q 2g _ n = r n _ 1 (g(«- 1 )/2 50 ) ; 

where m k = m k (c;logq) and T n (a,z) is in (|4.1.9p . By the formula of P 2g - n {m 2g -n)~ x Q 2g - n 
in the proof of Lemma 12.31 the latter equality means 

[cg{qW 2 ,z) ■■■ c-g+niqW^z)] [M 2g _ nA M 2g _^ 2 \ 

[sgiqi^y^z) ■■■ s-g+niqW^z)] [M 2g . nfi M 2g _ nA ] 

= [ S g(q("-l)/ 2 ,z) -.. S _ g+re _ 1 (Q(«- 1 )/2 )2 )]. 
These equalities follow from elementary identities 

c g _ {h+l) {q^\z) = c.g +n+h (q^/ 2 ,z), 

s g . {h+1) (q^- 1)/2 ,z) = -s_ g+n+h {q^/ 2 ,z) 
and the definition of M k j (1 ^ j < 4) in the proof of Lemma 12.31 □ 
Proof of (5). This is obvious by definition (|4.1.8|) . □ 
Proof of (6). Suppose that g^" 1 )/ 2 $C a < q n . By obvious equalities 

d d 
a—c k (a, z) = (-iz)s k (a, z), a—s k (a, z) = (-iz)c k (a, z), 

we have 
d 



a 



da 



x 



A q (a,z) 


z 


'-i " 




B q (a,z)_ 


~ 2 


-1 




g (a,z) 6 


s-i (a, 


z) ... 


















c g (a,z) c g -i(a,z) 



-g+n 



(a,z) 



v g (n). 



By definition (|2.1.2p of v g (n), we have P 2g _ n {m 2g _ n )v g {n) = Q 2g - n v g (n— 1). By looking 
at the (2g — n)-th row from the bottom, we have 

v gi n )[j\ = mig-n ■ v g (n)[j + 2g + 1 - n] for 1 < j < 2g - n. 
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This equality implies that 



s g (a,z) s g -.i(a,z) ■■■ s- g+n (a,z) 
••• c g (a,z) c g -i(a,z) 



c- g+n {a,z) 



v g (n) 



1 

1 



m 



2g-n 







m 2 g- n 

c g (a,z) c g -i(a,z) ■■ 


By using the identity 



c- g+n (a, z) 

s g (a,z) s g -i(a,z) 



-g+n 



(a,z) 



v g (n). 



-i 
-1 



we obtain 



1 

1 



d_ 

da 



in 







2g~n 
m 2 g-n 





"0 -l" 






1 





111 



2g-n 
m 2 g- 



1 

-i 



A q (a,z) 




"0 -l" 




B q (a,z)_ 


= z 


1 





111 



2g-n 







m 2 g- n 



A q (a,z) 
B q (a,z) 



for q^- 1 )/ 2 <^a< q n l 2 . This implies (6) by (HXTjl . 

Proof of (7). Note that ^(0) = A q (Q), since B q (z) is odd. We have 



□ 



lim 

o^-(qs)- 



A q (a,z) 
B q (a,z) 



1 

-i 

1 ' 

-i 



Cg(<l 9 ,z) 







2 




s g (q9,z) 



%(2<?)[1] 




by definition (I4.1.8p . Therefore, it is sufficient to show that 

v ff (2 5 )[l] =2]Tc fc + c 9 = A,(0). 

fc=0 

In order to prove this, we put 

S n = Po 1 Q Pi(m 1 )- 1 Q 1 ---P n (m n )- l Q n (n = 0, 1, 2, • • ■ ), 

where we understand rrik as parameters. The the size of S n is 2 x (2n + 4) by definitions 
of Pk(mk) and Qfc- We have 



So — Pq 1 Qo 
and find that the first row of S n has the form 



110 
1-1 



11 ••• 1 00 ••• 

n+2 n+2 



by the induction using the formula of -Pfc(m/%) 1 Qf c in the proof of Lemma 12.31 By 
definition (|2.1.2p . we have u s (2g)[l] = (the first row of S 2g -i) ■ v g (0). Hence we obtain 



v g (2g)[l) = 0_0_^_0) • v g (0) = 2Y J Ck + c g = A g (0). 

2g+l 2g+l k=0 



We complete the proof. 



□ 
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6. Proof of Theorem 12.41 
6.1. Preparations. We prepare two propositions for the proof of Theorem 12.41 

Proposition 6.1 Let 1 ^ a\ < clq ^ g 9 . 

(1) Assume that m q (a) ^ and m q (a)~ l ^ for every 1 ^ a ^ a^. Then there 
exists a 2 x 2 matrix valued function M(a±,aQ; z) such that all entries are entire 
functions of z, satisfies 



A q (ai,z) 
B q (a 1 ,z) 



M(ai,a ;z) 



A q (a ,z) 
B q (a ,z) 



(6.1.1) 



and det M(a±, ao; z) = 1. 
(2) Assume that m q (a) ^ and m^a)" 1 ^ for every 1 ^ a < ao- Then the matrix 
valued function M(ai,a; z) of (1) is left- continuous as a function of a and 



Ag(ai,z) 
B g (ai,z) 



A q (a,z) 
B q (a, z) 



lim M(ai,a; z) lim 

a— >o,q a^a,Q 
holds as a vector valued function of z € C. 
Proof. (1) Put A(a,z) = A q (a, z), B(a,z) = B q (a, z), m(a) = m q (a), and 

Then the system of (|4,l.lip in Theorem 14.21 is written as 



(6.1.2) 



-m{a) 




— 

da 



A(a,z) 
B(a,z) 



zJ{a) 



A(a,z) 
B(a,z) 



(1 ^ a < ao, z G C). 



By the assumption, m(a) and m(a) 1 are integrable on [ai,Oo]. Hence, we have 



A(ai,z) 
B(ai,z) 

A(a ,z) 
B(a ,z) 



A(a ,z) 
B(a , z) 



+ Z 



Ol 



A(h,z) 
B(t u z) 



A(a ,z) 
B(a ,z) 



dti 

ao rao 



oi Jti 



r° t / x^i 

U + W J( tl )-± + z 2 / J^JJ^J-i-i 

V Ja! H Jai Jti *2 *1 



J(h)J(t 2 ) 



A(t 2 ,z) 
B(t 2 ,z) 



dt 2 dt\ 



+z J 



0.0 rao rao 



at 



J {tl )J(t 2 )J(t 3 )^ d ± + ... 

13 12 El 



A(a ,z) 
,6(00,2;) 



where / is the 2x2 identity matrix. On the other hand, we have 



(6.1.3) 



J(h)-J(t k ) = {-!)* 





m(t2)rn(t4)-- ■m(f fc _ 1 ) 
. m(t 1 )m(t- A )---m.(t k ) 
■m(fi)m(t 3 )---m(f fc _i) 





m(fi)m(t3)---m(f fc ) 
m(t2)m(t4)---m(t fc _ 1 ; 




m(t 2 )m(t4)— m(t fc ) 
m(ti)m(t3)---»?i(t fe _ 1 ) 



if & = 2fc' + 1, 
if k = 2k'. 



Therefore, by taking 

C(a ,ai) 
and by using the formula 



sup{?n(a), m(a) ; a G [oi,ao]} 
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we obtain 



"0 



a ° , s , N dt k dt 2 dti 

J(ti)-J(tk)-r----r-r 



tfe-i 



fc! 



<^C7(ao,ai)*(log^ 



Ol 



for every 1 ^ i,j ^ 2, where [M]jj means the (i, j)-entry of a matrix M. This estimate 
implies that the right-hand side of (|6.1.3p converges absolutely and uniformly if z lie in 
a bounded region. 

Suppose that m(a) = m/ for a\ a ^ clq. Then 



I+z 



J(ti)-^+z 2 



ti 



J(t 1 )J(t 2 )^^+, 3 

12 *1 



ti ^2 



M J(tl)J(t2)J(t3) *** + . 



is equal to 

cos(z log(ao/ai)) — msin(z log(ao/ai)) 
^sin(zlog(a /ai)) cos(z log(a /ai)) 
and hence ()6.1.ip holds by taking this matrix as M(a\, oq; z). Therefore, if we suppose 
that m(a) = rrij / on [tj+\,tj) for a partition [oi,ao] = [ai,tfc_i) U • • • U [ii,ao] with 
*0 = a an d tk = 0-1, then we have (|6.1.ip by taking 

M(aQ,a 1 \z) 

cos(zlog(t fc _i/ai)) -m k sm(zlog(t k -i/ai)) 
^sin(zlog(t fe _i/ai)) cos(zlog(i fc _i/ai)) 

cos(zlog(t fc _ 2 /tfc_i)) -m fc _isin(zlog(t fc _ 2 /tfc_i)) 
^-^sm(zlog(tk-2/tk-i)) cos(^log(t fe _ 2 /tfc-i)) 



*3 *2 *1 



cos(zlog(ti/t 2 )) 
^sin^log^/^)) 

Moreover, det M(ai, ao; 2) 



-m 2 sin(zlog(ii/t 2 )) 
cos(> log(ti/t 2 )) 



cos(zlog(a /*i)) 
_^ T sin(zlog(a /ti)) 

1 is obvious by this definition. Now we complete the proof 



misinO log(a /*i)) 
cos(zlog(a Ai)) 



since m(a) is a constant on [q( n 1 ^ 2 , g™/ 2 ) for every 1 ^ n ^ 2g by definition (|4,1.6p . 
(2) By the above definition, M(ai,a; z) is left-continuous with respect to a, since m q {a) 
is left-continuous by definition (|4.1.6p . Because A q (a,z) and B q (a,z) are left-continuous 
with respect to a by Theorem 14.21 (4), we obtain (|6.1.2p from (|6. 1 . 1 j) . □ 

Corollary 6.2 Let m q (a) be of (|4.1.6f> . Assume that m q (a) ^ and m q (a)~ 1 ^ /< 



or 



every 1 ^ a < q 9 . Then we have 

co$(zlog(q n ' 2 / a)) 



A q (a,z) 
B q (a,z) 



E q (0) 



m 2g-n Sm 



(zlog(^/ 2 /a)) 



2g-n 

fc=i 



cos 



??? 



flog?) 

2 9 -(n+fe) sin (| lo g°) 



-m 29 „ n sin(zlog(g" /2 /a)) 
cos(zlog(g n//2 /a)) 

-m 29 „ (n+fc) sin(|logg) 
cos (f log g) 





"1" 








/or g^- 1 )/ 2 < a < a-™/ 2 . 



Proof. By Theorem 14.21 we can apply the proof of Proposition 16.11 (1) to a\ = a 



(q^- 1 )/ 2 <: a < 
tain the formula 



/2 ), a 
Proposition 6.3 Define 



q 9 and tk = q 



(g-k)/2 



(1 ^ k ^ g — n) and then we ob- 

□ 



Then we have 



K(a\ z, w) :- 



K(a; z, w) 



Eg(a, w)Eg(a, z) - E\{a, w)E\{a, z) 
2ni(w — z) 



A q (a, w)B q (a, z) - B q (a, w)A q (a, z) 



tt(z 



w 



(6.1.4) 



(6.1.5) 



ON ZEROS OF SELF-RECIPROCAL POLYNOMIALS 



21 



Moreover, if m q (a) and m q {a) 1 are integrable on [ai,ao], then we have 
K(a\; z, w) — K(ao; z, w) 



A q (a,w)A q (a,z) 



1 



m q {a) a 



da 1 
+ - 



7T 



ao 



B q (a, w)B q (a, z) m q {a 



«i 



for every z,w G 



da 



(6.1.6) 



Proof. We obtain (|6.1.5|) easily by substituting (|4.1.3p and into (|6.1.4|) . By the 

integration by parts together with (|4.1.1ip . we obtain 



ao ^T7 r < / v 1 da 

A(a, u;)/L(a, z) — — — 
ai m(a) a 



da 



A(a, w)B(a, z) 



<M) 



+ W 



ai 



B(a, w)B(a, z) m{a 



B(a,w)A(a,z) + w A(a,w)A(a, z) 



"i 

ao 



B(a, w)B(a, z) m(a 



1 



da 



ai 



da 
m(a) a 



Moving the second terms of the right-hand sides of the two equations to the left-hand 
sides, then adding both sides of the resulting two equations, and finally dividing both 
sides by (z — w), 



ao 



A(a, w)A(a, z) 



1 da 



ai 



m(a) a 



+ 



«0 . 



B(a, w)B(a, z) m{a) 



da 



{-A{a, w)B(a, z) + B(a, w)A(a, z)) 



w 



■rr(^K(ai; z, w) — K(ao; z, w)j . 



This implies (|6.1.6p . 



□ 



6.2. Proof of necessity. By Lemma |4~T| we see that all zeros of P g (x) lie on T and 
simple if and only if E q (z) is a function of class HB . Thus P g {x) has a multiple zero or a 
zero outside T if and only if E q (z) is not a function of class HB. Hence, for the necessity, 
it is sufficient to prove that E q (z) is not a function of class HB if ra2 9 _ n (c ; log q) ^ or 
m 2g-n(c j log <?)~ 1 ^ for some 1 ^ n ^ 2g. We will prove it in three steps as follows, 
there we put m 29 _ n = m 2g - n (c ; log q). 

Step 1. Firstly, we note that m 2g ^i > by (|5.0.2p . Therefore, if m 2g - n {c ; log q) ^ 
or m 2g ^ n {c ; log g) _1 for some 1 n ^ 2g, then there exists 1 ^ uq ^ 2g — 1 such 
that m 2g - n > and m^_ n > for every 1 ^ n < n , rn 2g _( no+1 )(c;\ogq) < or 
m 2#-(no+i)0s i 1°S <7) _1 ^ 0) A q (a,z) and B q (a,z) are defined for 1 ^ a < q n °l 2 and 



A q (z) 




~A q (l,z) 






B q (l,z) 



M(l,a;z) 



A q (a,z) 
B q (a,z) 



(6.2.1) 



holds for 1 < a < q n °/ 2 by applying (|6.1.ip to (a\,ao) = (l,o). 

Step 2. Let no be the number of Step 1. Suppose that E q {a, z) = A q (a, z) — iB q (a, z) 
is not a function of HB for some 1 < oq ^ q n °/ 2 , that is, E q (ao,z) has a real zero for 
some 1 < ao ^ q n °/ 2 or | .©5(010, 2) | ^ \E q (ao,z)\ for some Imz > and 1 < ao ^ q n °^ 2 . 
Here, we understand A q (ao,z) and B q (ao,z) in the sense of left-sided limit a — > a$ if 
a = «f°/ 2 (see (J6JL2])). 

If E q (a, z) has a real zero for some 1 < ao ^ q n °^ 2 , then A g (ao, 2) and B q (ao, z) have 
a common real zero, since they are real valued on the real line. Therefore, (|6.2.ip and 
det M(l, ao; z) = 1 imply that A q (z) and B q {z) have a common real zero. Hence E q {z) 
has a real zero and thus E q (z) is not a function of class HB. 

On the other hand, we assume that E q (a, z) has no real zeros for every 1 < a ^ q n o/ 2 
but it has a zero in the upper half plane for some 1 < ao ^ qW 2 . By (|4,1.8p and 
Theorem 14.21 (3). E q (a,z) is a continuous function of (a, z) G [l,q™ () / 2 ] x C. Therefore, 
any zero locus of E q (a,z) is a continuous curve in C parametrized by a G [l,g n0//2 ]. 
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Denote by z a C C a zero locus through a zero of E q (ao,z) in the upper-half plane, 
that is, E q (a,z a ) = for every 1 ^ a ^ q n °l 2 . If Im(z ai ) < for some 1 ^ a± < ao, 
then Im(,z a2 ) = for some a\ < a 2 < ao- This implies that E q (a 2 ,z) has a real zero 
at z = z a2 . This is a contradiction. Therefore, Im(z a ) ^ for every 1 ^ a < ao, in 
particular Im(zi) ^ 0. This implies E q (z) = E q (l,z) is not a function of class HB. 

If E q (a, z) ^ for Imz > but \E q (a , zq)\ ^ |-E ? (aoi z o)| for some 1 < a ^ g n o/ 2 and 
Im(zo) > 0, then it derives a contradiction. Because A q (a, z) and B q (a,z) are bounded 
on the real line by definition (|4.1.8p . E q (a,z) is a function of the Cartwright class |141 
the first page of Chapter II]. Therefore, we have the factorization 



E, 



q{ao,z) 



C lim 



n 

\p\<R 



1 



(see |14| Remark 2 of Lecture 17.2]). Here Im(p) < for every zero of E q (aQ,z) by the 
assumption. Hence, we have 

Eg(ao,z) 



E q (a ,z) 



lim 



n 

\p\<R 



< 1 for Imz > 0. 



This contradict to the assumption \E q (ao, zq)\ ^ \E q (aQ, zo)\. 

Step 3. For the number uq of Step 1, one of the following case occurs: 



(i) m 2 g- 



or m n 



-v'io+i) ' 
(ii) rn 2g _ {no+1) < and m£_ {nQ+1) + 0. 

We prove that E q {z) is not a function of HB even if whichever occurs. Considering 
the argument in Step 2, we assume that E q (a,z) is a function of class HB for every 
1 < a ^ q n o/ 2 in both cases. 

Case (i). Suppose that m2 9 _( no +i) = 0. Then, by definition (|2,2.2p . we have 

v g (n )[l] + v g (n )[2g - n + 1] = 0. 
This implies that the function 

A q (a,z) = l -(T n ,(q n °l\z)v g (n ))[l) ((^/^-o^ 



lim 

a->(g n o/ 2 )- 



is of exponential type whose mean type is at most (g — 1 — (no/2)) log q. If 

^(n )[2o - n + 2] + v g (n )[4g - 2n + 2] = 
for the denominator of (|2.2.2p . then the function 



lim B q (a,z) 

a^(q n a/ 2 )- 



;( T n {q 



n / 2 



Z V, 



("o))[2] 



is also of exponential type whose mean type is at most (g 
other hand, we have 



1 - (n /2))logg. On the 



A q (l,z) 
B q (l,z) 



lim 

a^(<j"0/ 2 )- 



M(l,a;z) 



<f>(zj 
Hz) 



by (|6.1.2p . Here, the entries of the left-hand side are entire functions of exponential 
type with mean type glogq, while entries of the right-hand side are entire functions of 
exponential type with mean type at most (g — 1) log q from the construction of M(l, a; z) 
in the proof of Proposition 16.11 (1). This is a contradiction. Hence, it must be 

v g (n )[2g - n + 2] + v g {n Q )[Ag - 2n + 2] ^ 

and ip(z) is a function of exponential type whose mean type is just (g — (no/2)) logg. 

By the assumption, E q (a, z) is a function of HB at a = q n °/ 2 (in the sense of left-sided 
limit). Therefore, (j>{z) and ip(z) have only real zeros and their zeros interlace. However, 
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by |14| Theorem 1 of Lecture 17.2], main terms of asymptotic formulas for the number 
of (real) zeros of (f>(z) and ip(z) in [— T, T] C M are strictly different. In particular, their 
(real) zeros can not interlace. This is a contradiction. Hence E q (a,z) is not a function 
of HB for some 1 < a ^ q n °/ 2 which implies that E q (z) is not a function of HB by Step 
2. The case of fn^g-no = is proved by a similar way. 

Case (ii). In this case, we can assume that E q (a,z) is a function of class HB for 
every 1 < a ^ q>( n o+i)/2^ s i nce we only used vri2g- n 7^ and fn^g-n 7^ in Step 2. Put 
0l = g n /2 ; aQ = ^(n +i)/2 _ q n /2y 2 and m 2g _ {no+1) = -m < 0. Then, for every 
cti ^ cl ^ ao, we have m q {a) = —m by (|4.1.6p and find that E q (a,z) generates the de 
Branges space B(E q (a, z)) which is the Hilbert space of all entire functions F{z) such 
that L \F(x)/E q (a, x)\ 2 dx < oo and F(z)/ E q (a, z), F(z)/E q (a, z) are functions of the 
Hardy space H 2 in the upper half-plane (see [5| §19] and |18[ Proposition 2.1]). 

We have K(ao; z, z) > K(a±; z, z) by applying (|6.1.6p to z = w with m q {a) = —m < 0. 
Therefore, it follows that for every / G B(E q (ai, z)) 

\f(z)\ 2 ^\\f\\ 2 ai K( ai ;z,z) <\\ff ai K(a ;z,z) 

by [H Theorem 20], where || • || ai is the norm of B(E q (a\, z)). By applying this to the 
function 

:= E q (a 1 ,z)-E q (a 1 ,iy ) 
z - «yo 

which is a function of B(E q (a\, z)) by Lemma 3.3 and Lemma 3.4 of [5], we obtain 

|2 \E q (a ,iy)\ 2 - \E\(a Q ,z)\ 2 2 \E q (a , iy)\ 2 



\g{iy)? < htM^vM = \\g\ V 9V *'\ ™ " < llfllSi- 



47ry 1 47ry 

By E q (a, z) = A q (a, z) — iB q (a, z) with (|4.1.8p . we see that 



n ° 2 1)v < \g{iy)\ < 2/ 1/2 |^(a ,Ji/)| «y V V 9 as y -> +oo. 



y q 



This is a contradiction. Hence E q (a, z) is not a function of HB for some 1 < a ^ ao and 
it implies that E q {a, z) is not a function of HB by Step 2. □ 

6.3. Proof of sufficiency . Suppose that m^g-nisk ; log g) > and m2g- n (c;logg) 1 > 
for every 1 ^ n ^ 2a. Then m q (a) and m g (a) _1 are both integrable on [1, q 9 ) and positive 
real valued. Therefore, by applying (|6.1.6p to (a±, Oo) = (a, b) and then by tending b to 
q 9 together with Theorem 14.21 (7) and (|6.1.4p . we have 

0<-/ \A q {t,z)\ 2 —r-- + - B q (t,z) 2 m q (t)- 

^ x u(h x l^(a,z)| 2 -|4(a,z)| 2 

= A (a; z, z) — hm A (o; z, z) = — - — - 

b^qS 47rlm z 

for every 1 ^ a < q 9 if Imz > 0. Thus E q (a,z) is a function of class HB for every 
1 ^ a < q 9 . In addition, we have 



lim M(a, 6; z 







B q (a,z) 

for every 1 ^ a < (7 51 . Here E q (0) = A q (0) and 

^ g (2g)[i] _ 4/(0) 



with lim det M(a, b;z) = 1 



7^ m (c; logo) 



w fl (2 5 )[2] w fl (2 5 )[2] 

by the proof of Theorem H3I (7). If 4,(0) = v 9 (2o)[l] = u ff (2g)[2] = 0, then A,(a,z) 
and B q (a,z) are both identically zero for q( 2 5 l_1 )/ 2 ^ a < q 9 by (|4.1.8p . It implies 
that E q (z) = A q (z) = B q (z) = by Proposition 16.11 together with the assumption. 
Therefore, A q (0) 7^ and it implies that A q (a, z) and B q (a,z) have no common zeros 
for every 1 ^ a < q 9 . Thus E q (a,z) has no real zeros. As a consequence E q (a,z) is a 
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function of class HB for every 1 ^ a < q 9 . Hence A q (z) has only simple real zeros by 
Eg (z) = E g (l,z) and Lemma Q (2). □ 

7. Proof of Theorem 12. 7L 12.91 14.41 and Corollary 14.51 

Firstly, we prove Theorem 12.71 by using Theorem 14.41 as well as the proof of Theorem 
12.41 However, we omit the proof of Theorem 14.41 because it is proved by a way similar to 
the proof of Theorem l4.2l in Section [5j Also, we omit the proof of Theorem I2.8| because it 
is just a rewriting of Theorem 12. 71 by (I2.3.3P and ()2.3.4p . Successively, we prove Corollary 
14.51 by using Theorem 12.71 and Theorem 14.41 as well as Corollary 14.31 Finally, we prove 
Theorem [ 



7.1. Proof of Theorem 12. 71 We start from the following lemma. 

Lemma 7.1 Let E q ^(z) be in (|4.2.ip . Then, all zeros of P g (x) lie on T if and only if 
E n ^(z) is a function of class HB for every uj > 0. 

Proof. By definition (|4.1.1|) . all zeros of P g {x) lie on T if and only if A q (z) has only real 
zeros. Suppose that A q (z) has only real zeros (allowing multiple zeros). Then E qiJ (z) 
satisfies inequality (|3.0,5p for every uj > 0, because we have 





2 


A q (z — iuS) 


2 

-n 

p 


(x 


- p) + i{y-uj) 


E q,u>{z) 




A q (z + iuj) 


(x 


- p) + i{y + uj) 



n 



4ujy 

1 < 1 



(x - p) 2 + (y + uj)'' 



for z = x + iy with y > by using the factorization (|4.1.5p . Moreover, E qu (z) has no 
real zeros for every uj > by definition ()4.2.ip and the assumption. Hence E q<u (z) is a 
function of class HB for every uj > 0. 

Conversely, suppose that E qw (z) is a function of class HB for every uj > 0. Then all 
zeros of A q ^{z) and B qtUJ (z) are real, simple, and they interlace by Proposition 13.31 In 
particular, A qtUJ (z) has only real zeros for every uj > 0. Hence A q {z) = lim w _;.o+ A q)W {z) 
has only real zeros by Hurwitz's theorem in complex analysis ( |16| Th. (1,5)]). □ 

Proof of necessity. By Lemma 17.11 P g (x) has a zero outside T if and only if E q ^ w {z) 
is not a function of class HB for some uj > 0. Hence it is sufficient to prove that if there 
exists ujq > such that m2 9 _ n (c; q 1 ^ ) ^ or m2 5 _„(c; q ul0 )~ 1 ^ for some 1 ^ n ^ 2g, 
then E qtUJo (z) is not a function of class HB. This is proved by a way similar to the proof 
of Section 16.21 by using Theorem 14.41 instead of Theorem 14.21 □ 

Proof of sufficiency. Let uj > 0. Suppose that m2 g - n {c ; (f) > and m,2 g - n (c ; q^)^ 1 > 

for every 1 ^ n ^ 2g. Then it is proved that E q ^{z) is a function of class HB by 
a way similar to the proof of Section 16.31 by using Theorem 14.41 instead of Theorem 
14.21 Therefore, by Lemma 17.11 if m2 9 _ n (c; g w ) > and m2 9 _ n (c ; q^)^ 1 > for every 

1 ^ n ^ 2g and uj > 0, then all zeros of P g (x) lie on T. □ 

7.2. Proof of Corollary 14.51 If all zeros of P g (x) lie on T, then A qtUJ (z) has only simple 
real zeros for every uj > by Proposition 13.31 and Lemma 17.11 In particular, E q ^(0) = 
Aq t u(0) 7^ 0, since A qui (z) is an even function having only simple real zeros. Hence, 
H q<u (a) defines a canonical system on [l,q 9 ) having the solution (A q ^(a, z)/E qtUJ (0), 
B q ^ U} (a,z)/E q ^ UJ (0)) by Theorem 12.71 and Theorem 14.41 Conversely, if H qoJ (a) defines a 
canonical system for every uj > 0, then m2 g - n (c.',Q UJ ) > and m2 9 _ n (c ; q u ')~ 1 > for 
every 1 ^ n ^ 2g and uj > by Definition 13.11 Hence, all zeros of P g (x) lie on T by 
Theorem 12.71 
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7.3. Proof of Theorem [2JH We define m 29 -„(c;<fO for 1 < n SC 2g by (|2~TTB and 

a <7,w(°) 
_u; _1 a gA ,(0) 

7^ _)_ 



(|2.1.2p starting from the initial vector v g (0) 
of (|2.3.ip . In addition, we take rh,2g{c;q u) ) := U) 



qguj _ q-gu 
rh2g-n(c;q u ) = uj ■ m 2g - n (c;q UJ ), 
for every 1 ^ n ^ 2g, where m2g- n {c]q w ) is of (|2.3.2p . This implies 

Rn(c;(f) = B^(c;q w ) 
for every 1 ^ n ^ 2g if we define i? n (c ; g w ) by 



where a goj (0) is the vector 
. Then we have 



Rn(c;q u 



( J 

n 

j=0 

j 



m 2g _(2j+i)(c;g a ') 
fri2g-2j(c;q u ) 



if n = 2J + 1 > 1, 



n m 2g -( 2j+2 )(c;<n ifn = 2J+2 ^ 2 

■ L ^"l23-(2i+l)(c;^) 



Therefore, we obtain 



qgu _)_ _ 
m2g-n{c;q") = — ::= ^Rn-i{c;q")R n {c;q u 



qguj _ q-g<+> 
1 / 1 



5 V wl °g9 

Hence, for Theorem 1 2. 9 1 it is sufficient to prove that 



lim Rn(c;(f) = R n (c) (l<n<2 5 ). 

g<"->-l+ 



By the definitions of ^j-nfc;^) an d TO2 9 _n(c;log(/), this equality follows from the 
following formula of v g (l): 



lim P 29 -l("l2g-l) 1 Q2g-l 



a 9,o;(0) 
~S,"( ). 



f2ff-l( m 2ff-l) 1 Q2g-l 
where the right-hand side is the vector of f|5.0.3j) . Put 

P2g-l(m2g- n {c; q° J ))~ 1 Q2g-l 



a s (0)" 




a 5 (l)" 


b 5 (0)_ 




b s (l)_ 



a g,o>(l) 
b ff ,a,(l). 



a 9,w(0) 

w _1 a 9ia ;(0) 
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Then, by using the formula of Pk(mk)~ 1 Qk hi the proof of Lemma |2,3| we have 

2 cosh (g log q LU )cQ 
(cosh((<? — 1) log q^) + oj m 2g -\ sinh((g — 1) log <p))c\ 
(cosh((# - 2) log q u ) + oj~ x m 2g -\ sinh((# - 2) log q u ))c 2 



(cosh (log q u 
(cosh (log q u 



oj 



OJ 



~ 1 m 2g -i sinh(log <f))cg-i 
' 1 m 2g -i sinh(log <f))c g -\ 



(cosh((<? — 2) log q w ) — oj 1 rri2 g -i sinh((g — 2) log q^))c 2 
.(cosh((# - 1) log q^) - uj- 1 m 2g -i sinh((# - 1) log q u ))ci. 

2 oj~ 1 sinh(g log q u )c$ 
m 2g l -i{ cos]n -((9 ~ l)logg") + uj- 1 m 2 g-ismh((g - l)log^))ci 



■2g 



, 1 _ 1 (cosh(( 5 - 2) log <F) + uj- x m 2 g-i sinh(( 5 - 2) log <f))c 2 



m 



2.y 



,_i(cosh(logg a) ) + oj 1 fn 2g -\ sinh(log (f))c, 



■9-1 



m 2g-l C 9 



m 



2.y 



x (cosh (log cf) - oj 1 m 2g -\ sinh(logg"))c 5 _i 



m 



2 H 



, 1 _ 1 (cosh(( 5 - 2) log q") - oj~ x m 2 g_ x sinh(( 5 - 2) log g w ))c 2 



m 2o-l( cosn ((s ~ l)logg") -w 1 m 2 3_i sinh((5 - l)log^))ci 



with 



By using 



m 2g -i = OJ 



j9u _|_ q-gu 



-gu 



a;coth(<? log q w ). 



lim 

lim 

i->0+ 



cosh((<7 — k)x) + coth(gx) sinh((g — k)x) 

cosh((<? — k)x) — coth(gx) sinh((g — k)x) 
logg 



2g-k 
9 



(O^k^g), 



9 



lim m 2 1 = lim tanh(^x) = glog<7 

a;->0+ y z->0+ X 



for x = log q u we obtain 



lim 







\0)' 


Pg,M). 




b 9 (l)_ 



On the other hand, it is easy to see that 



lim m 2g (c;0 = — : 



"J2 S (c, logg). 



Hence lim (? ^^ 1 + R n {c'-,q UJ ) = R n (c) f° r every 1 ^ n ^ 2g. Because R n {c;q u ) are rational 
function of g", we obtain the second formula of Theorem 12.91 Therefore, we obtain 

m 29 _ n (c;g w ) = m 2g (c; q UJ )R n ^ 1 (c; q UJ )R n (c; q") 
1 



ug \ogq 

This implies the third formula of Theorem 12.9 



+ 0(log<n i4-i(c) + 0(log^) i? re (c) + 0(logg 



□ 
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7.4. Remark on Theorem 12.91 We have 

E q ^{z) = A q (z) - iujB q (z) + O z (uj 2 ), 



, w (z) = A q {z) + O z (uj z ), B q ^(z) = ujB q {z) + O z {uj 6 ) 

as uj — > + if z lie in a compact subset of C. Therefore, it seems that E q ^{z) is similar 
to E q {z) = A q (z) — iB q (z) for small uj > 0, but there is an obvious gap after taking the 
limit uj — > + . To resolve this gap, we consider 



E q ,u{z) := A qiU (z) B qtU] {z) 



UJ 



Then, we have 



E q M = A q {z) -iB q {z) + O z {uj 2 ) = E q (z) + O z (l 



A 



(z) := \{E q ,u{z) + E\jz)) = A q>UJ (z) = A q {z) 



+ O z (uj 2 ), 



) - ElJz)) = ±B qjU (z) = B q {z) + O z (uj 2 ) 



as uj — > + if z lie in a compact subset in C. Hence E q ^(z) "recovers" E q {z) by taking the 
limit uj —7- + . The initial vector v g (0) of Section \7. 31 is chosen so that it corresponds to 
A qtU (z) and B qtLJ (z). This is a reason why Theorem 12 . 9 1 holds . In spite of this advantage, 
we chose E qjUJ (z) not E q ^(z) to state results in Section [2.31 One of the reason is the 
simple formula E qtU} (z) = A q (z + iuo). Comparing this, E qtU) (z) has a slight complicated 
formula 

~ , , _ A q (z + ioj) + A q (z — ioj) A q (z + iuj) - A q (z - iuj) 
E q ,M) ~ 2 + 2uj • 

We do not know whether an analogue of Lemma [7.11 holds for E q ^ w {z). However, if such 
analogue holds, we may obtain results for E q ^(z) analogous to results in Section [2.31 

8. Concluding Remarks 

By Corollary 16.21 and the proof of Section 16.31 a self- reciprocal polynomial P g {x) of 
(jl.O.ip has only simple zeros on the unit circle if and only if there exists 2g positive real 
numbers mo, ■ ■ ■ , fnig—\ such that 



n=l 



(x + 1) 

-1 



-im 2g _ n {x 



im 2g -n{x - 1) 
1) (x + 1) 





Y 








since P g (l) = E q (0) = A q (0). Compare this with the factorization 



2XjX + 1) (Xj G 



P 9 (X) = P g (0)H(x 2 

3=1 

As in Section [21 we have a simple algorithm to calculate uiq, • • • , m,2 9 _i from coefficients 
Co, • • • , c g , but it is not for Ai, • • ■ , X g , since the Galois groups of a general self- reciprocal 
polynomial P g (x) is isomorphic to <5 g x (Z/2Z) 9 . In addition, it is understood that 
the positivity of uiq, ■ ■ ■ ,l7i2g—i is equivalent to the positivity of a Hamiltonian, but a 
plausible meaning of |Aj| < 1 and Aj 7^ Xj (i 7^ j) is not clear. 

Before concluding this paper, we remark two important remaining problems. 

(1) As mentioned in the end of Section [2.2[ general closed formula of m2 S „„(c, log q) 
or R n (c) is not yet obtained even for conjectural one. Therefore, the discovery of such 
formula is desirable. In fact, we need a simple closed formula of m2 9 - n (c, logg) or R n (c) 
for the convenience of actual applications of Theorem 12.41 or 12.61 The following formula 
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of R n (c) may suggests the existence of certain "nice" general formula of R n (c), although 
an expected result is not yet clear. We have 

9-1 9 

p g( x ) = Yl c k( x29 ~ k + xk ) + c 9 x9 = c o I]V " 2X i x + 1 )' 

k=0 j=l 

where 2\j (1 ^ j ^ g) are zeros of the Chebyshev transform of P g (x) ([S])- By using 
Ai, • • • , X g , we obtain the following formula of R n (c): 

•5 = 1) J? 2 (c ,ci) 



l-\t 



l + Xi 

• g = 2, R n = R n {c , a, c 2 ) (2 < n < 4), 

_ (i-Ai) + (l-A 2 ) _ o (l-Af) + (l-Aj) _ (l-A 1 )(l-A 2 ) 

^ 2 (l + A 1 ) + (l + A 2 )' K3 (Ax-A 2 ) 2 ' U4 (i + Al )(l + A 2 )' 

• 5 = 3, Rn = R n (c ,ci,c 2 ,c 3 ) (2 < n ^ 6), 

(1-A 1 ) + (1-A 2 ) + (1-A 3 ) q (1-Af) + (1-A|) + (1-A|) 

-fi2 — ~ — — — — r^r, ti,z 



(1 + Ai) + (1 + A 2 ) + (1 + A 3 )' ° (Ai-A 2 ) 2 + (A 1 -A 3 ) 2 + (A 2 -A 3 ) 2 ' 
„ Ei^< K 3a-^)(l-A J )(A,-A J ) 2 El ^ <K3 (l-A 2 )(l-A 2 )(A J -A J )^ 



(1-Ai)(l-A 2 )(l-A 3 ) 
(1 + Ai)(l + A 2 )(l + A 3 )' 

These formulas look simpler than formulas of R n (c) in the end of Section f2.2l 



R 



Remark After completing the first version of this paper, simple closed formulas of 
m2g-n(c, logg), m 29 _ n (c, q u ), R n (c) and R n {c]q w ) are obtained. See |22j for details. 

(2) There are several important classes of self-reciprocal polynomials with real coeffi- 
cients. Here we mention two of them. The first one is zeta functions of smooth projective 
curves C/¥ q of genus g: Zc(T) = Qc(T)/((l — T)(l — qT)), where Qc(T) is a poly- 
nomial of degree 2g satisfying the functional equation Q C {T) = (q 1 / 2 T) 2 9Q c {l/(qT)). 
Hence Pc(x) = Qc{l~ 1 ^ x ) is a self-reciprocal polynomial of degree 2g with real coeffi- 
cients. Weil |23| proved that all zeros of Pc(x) lie on T as a consequence of Castelnuovo's 
positivity for divisor classes on C x C . The second one is polynomials Pa{x) attached 
to n x n real symmetric matrices A = (a^.j) with |ajj| ^ 1 for every 1 ^ K j ^ n 
(no condition on the diagonal): Pa(x) = S/uJ={i 2 ■■■ n} Tliel j&J a M> wnere I U J 
means a disjoint union. Polynomials Pa(%) are obtained as the partition function of a 
ferromagnetic Ising model and they are self-reciprocal polynomials of degree n with real 
coefficients. The fact that all zeros of any Pa{x) lie on T is known as the Lee- Yang circle 
theorem |12j . Ruelle |19| extended this result and characterized polynomials Pa(x) in 
terms of multi-affine polynomials being symmetric under certain involution on the space 
of multi-affine polynomials |20] , 

It seems that a discovery of arithmetical, geometrical or physical interpretation of the 
positivity of m 25 _ n (c, ; log q) or R n (c) (for some restricted class of polynomials) is quite 
interesting and important problem. Such philosophical interpretation may contribute to 
find a simple closed formula of m 29 _ n (c, ;logg) or R n (c). 
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